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Abstract: A finite set statistics (FISST)-based method is proposed for multi-target 
tracking in the image plane of optical sensors. The method involves using signal amplitude 
information in probability hypothesis density (PHD) filter which is derived from FISST to 
improve multi-target tracking performance. The amplitude of signals generated by the 
optical sensor is modeled first, from which the amplitude likelihood ratio between target and 
clutter is derived. An alternative approach is adopted for the situations where the signal noise 
ratio (SNR) of target is unknown. Then the PHD recursion equations incorporated with 
signal information are derived and the Gaussian mixture (GM) implementation of this filter 
is given. Simulation results demonstrate that the proposed method achieves significantly 
better performance than the generic PHD filter. Moreover, our method has much lower 
computational complexity in the scenario with high SNR and dense clutter. 

Keywords: finite set statistics; signal amplitude information; PHD filter; multi-target 
tracking; Gaussian mixture 



1. Introduction 

Optical sensors have been widely applied both in important military and civil areas due to their 
properties of long detection range, high concealment ability and large coverage area. Since there is 
usually a long distance between targets and sensor which generates a low SNR and dense clutter 
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scenario, multi-target tracking in the image plane of optical sensors is a very difficult problem. In 
multi-target tracking, the aim is to estimate the number of a set of targets and the state of each target 
from a set of measurements received. However, due to the variation of targets number with time in the 
field of view of the sensors and the existence of miss detection and dense clutters, multi-target tracking 
in image plane remains a challenging problem. 

In addition to the location measurement, amplitude information has been proven to improve tracking 
performance [1-5]. One of the pioneering techniques was the approach proposed by Colegrove, Lerro 
and Bar-Shalom [1,2] where the probability data association (PDA) filter utilizing target amplitude 
was applied in the context of single-target tracking. The target amplitude has also been incorporated in 
the multiple hypothesis tracking (MHT) framework [3] and Viterbi data association scheme [4]. More 
recently, the significance of target amplitude has been explored for data association of closely spaced 
targets in [5]. Although significant progress has been made recently, the approaches mentioned above are 
based on data association technique, which requires expensive computational cost in most circumstances. 
Therefore, the traditional data association approach is not the optimal option for multi-target tracking 
in image plane of optical sensor in the scenarios with time-varying targets number, low SNR and 
dense clutters. 

A good candidate is the emerging Bayesian approach in the framework of finite set statistics 
(FISST) proposed by Mahler [6]. FISST provides a set of mathematical tools that allows direct 
application of Bayesian inferencing to multi-target problems. PHD proposed by Mahler [7] is a 
computation tractable approximation for the optimal multi-object Bayes filter based on Finite Set 
Statistics (FISST). Operating on the single-target state space, the PHD filter avoids the combinatorial 
problem that arises from data association and thus leads to superior performance in comparison with 
traditional MHT algorithm in multi-target tracking [8,9]. These features render the PHD filter 
extremely attractive to many researchers. The primary PHD applications can be found in [10-13]. 
In [14], Clark et al. incorporated target amplitude into a PHD filter for the first time in multi-target 
tracking, which improves the tracking performance. However, their amplitude is modeled for radar and 
sonar application and the computational complexity analysis of their method is absent. 

In this paper, we propose a FISST based method which uses signal amplitude information in a PHD 
filter for multi-target tracking in the image plane of optical sensors. Based on analyzing the imaging 
characteristics of the optical sensor, we model the signal amplitude and incorporate it into the PHD in 
the form of amplitude likelihood ratio. In the situation where the SNR of target is unknown, we then 
present an alternative method based on this amplitude model. Simulation results demonstrate a 
significant improvement of the proposed method in tracking performance. Furthermore, the 
computational complexity of our method in the scenarios with different clutter density and SNR is 
also discussed. 

This paper is organized as follows: Section 2 models the signal amplitude generated by the optical 
sensor. Section 3 demonstrates how to incorporate the amplitude information (for known and unknown 
SNR case) into a PHD filter and then the GM implementation for this filter is given. Section 4 presents 
the simulation results that validate the proposed method. The conclusions are given in Section 5. 



Sensors 2012, 12 



2922 



2. Amplitude Measurement Model 

2.1. Amplitude Likelihood Ratio 

In an optical multi-target tracking system, the original images taken by the sensor are usually 
processed by background suppression before being further used for target tracking. Assuming the noise 
is addictive, we define the SNR d of the residual image as [15-17]: 



d = s / 



(1) 



where s is the mean signal value of target, a is the standard deviation of the residual image. Each 
measurement from the image consists of the two-dimensional position vector z in the image plane and 
the corresponding amplitude a > 0, that is, a measurement vector has the form z:=(z r ,a) r . For 
simplicity, we assume that the value of the target signal has no spreading in the image plane. We 
assume that the noise is Gaussian, then the probability densities of the amplitude of the false alarms 
and the target po(a) and p\(a\d) can be written as [18]: 
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This leads to the probabilities of false alarm and detection p FA and p T D (d) with a detection 
threshold x, 

Pfa = J" Po ( a )da = erf(r/(7) (4) 

/• CO 

p T D (d) = J Pi(a\ d)da = erf(r/cJ - d) (5) 
dy is the probability error function. According to Equations (4) and (5), 



where erf(pt) = -=L= f + °°exp 
we have: 
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e Tf- l (p* FA )- e rf- l [p* D (d)] = d 



(6) 



In target tracking application we usually choose p r FA as a fixed value. Given p FA , the threshold x 
can be calculated via the inverse form of Equation (4) with the parameter a, and the probability of 
detection p T D (d) can be calculated via Equations (5) or (6) for the target with SNR = d. Table 1 lists 
the values of p T D (d) for different SNR values d and specified p FA with a setting to be a normalized 
value, i.e.,o= 1. 
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Table 1. p T D (d) under different SNR and p T combinations. 



d 



P F A 
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5 x 1(T 5 


0.5436 


0.6864 


0.9825 


0.9991 


1.0000 


1 x 1(T 4 


0.6106 


0.8999 


0.9887 


0.9995 


1.0000 


5 x 10~ 4 


0.7610 


0.9563 


0.9966 


0.9999 


1.0000 


1 x 1(T 3 


0.8185 


0.9719 


0.9982 


1.0000 


1.0000 



When the target SNR d is known, we can get our amplitude likelihood functions for the false alarm 
and the target as: 

c T a{a) = -^—p 0 (a), a>T (7) 
Pfa 

sl{a\d) = —^—p l (a\d), a>T (8) 
p D (d) 

The amplitude likelihood ratio given a threshold x is defined as [19]: 

c fl (o) 

We use the notation c a (a) and p a {a\d) to denote the case where x = 0, then we have: 

c a (a) = p 0 (a) (10) 

Pa{a\d)= — (11) 
p Q (a) v ' 

From Equations (5) and (9) we can see that the calculations of p T D (d) and p x a (a\d) rely on a 
specified known target SNR, however, this requirement cannot be satisfied in most practical tracking 
systems. We adopt an alternative approach to circumvent this issue next. 

2.2. Method for Unknown SNR 

When the SNR of target is unknown, one straightforward approach would be to estimate the 
unknown parameter d from the measurement amplitudes a. However, this approach requires a large 
number of measurements from the target to achieve an accurate estimate of d. Furthermore, due to the 
unknown association between measurements and targets in multi-target environment with clutter, the 
approaches of estimating d usually fail. Similar to the idea introduced in [14], we adopt an alternative 
approach where we do not attempt to estimate d at all. Instead we marginalize out the parameter d over 
the range of possible values and find a probability of detection for p T D and a likelihood ratio for p\ 
that is not conditional on d. 

Consider that we always have some prior information about the targets being tracked, so we assume 
that p(d), defined on the possible SNR values [d\,di\, gives the expected probability distribution of 
SNR values. Since the amplitude distribution in Equations (2) and (3) are symmetric thus have no 
biases in high or low SNR targets, the reasonable choice for p(d) is the uniform distribution U\d\,dj\. 
Then we define the probability of detection and amplitude likelihood ratio where SNR is unknown as: 
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p T D=l d 1 p(y)p T D(y)dv (12) 

P T a{a) = \ 2 P(v)pl{a\v)dv (13) 



From Equations (5) and (12) we have: 



4-4* - < 14 » 

Note that the p T D over the marginalized region [d\, fife] can be computed with numerical integration 
offline since it does not need to be computed at each iteration. The computation of p T a (a) as in 
Equation (13) can be simplified and will be presented in Section 3.2. 

3. PHD Filter with Signal Amplitude Information 

Suppose that at time k there are N k target states x k l ,---x k Ni , each taking values in a state space 
X cR"' ; and M k measurements (detections) z kil ,---z kM , each taking values in the observation space 
Z c R" . In PHD filter, a multi-object state and a multi-object observation are represented by RFS: 

x k ={x k ^-x ktNk ]e¥ (X) 

Z k={zk,l>---Zk,M k } G¥ ( Z ) 

where F(X) and F(Z) are the finite subsets of X and Z, respectively. The state x = (x,y,x,y) T of each 
target contains the position (x,y) T and velocity (x,y) r in the image plane, while the measurement zis 
defined in Section 2. We assume that each target follows a linear Gaussian dynamical model and the 
sensor has a linear Gaussian measurement model, i.e., 

fk\k-\( x I *') = N {x;F k _ x x,Q k _ x ) (16) 

L z (z\x) = N (z;H k m,R k ) (17) 

where N (■;/«, -P) denotes a Gaussian density with mean m and covariance P, F k .\ is the state transition 
matrix, Q k .\ is the process noise covariance, H k is the observation matrix, and R k is the observation 
noise covariance. 

3.2. The PHD Recursion with Amplitude Information 

We abbreviate the PHD filter incorporated with amplitude information as AI-PHD filter. Next we 
derive the prediction and update equations of AI-PHD filter based on the amplitude likelihood ratio 
given by Section 2. For simplicity, we do not consider target spawning in this paper. 

Step 1. Prediction: The prediction equation of AI-PHD filter is the same as generic PHD filter since 
their state vector and state transition matrix are the same, i.e., 

D k\k-\ (*) = 7k (*) + f (Ps,k (x')f k \ k -i (x | x))D k _ nk _ x (x')dx' (18) 
where y k (x) is the birth term for new targets, p Sk (x) is the probability of target survival, f k ^ k _ l (x\x) is 
the transition density and D k _ Vik _ l (x') is the PHD at time k— 1. 
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Step 2. Update: The update equation is changed when incorporated with amplitude information. 
Analogized to the update equation of generic PHD filter in [7], we have the update equation of our 
AI-PHD filter as 

D k\k{x) = L 2i (x)-D k]k _ x (x) (19) 
where L 2 (x) is the pseudo-likelihood function as 

% t (*) = 1 - Id (*) + Id (*) £ ; v /-x * r Fl (20) 

D k\k-\ [Pd L z1 = j{Pn { x ) L i {x)D k{k _ l (x)}dx (21) 

where X and V are the clutter density and area of image plane of optical sensor respectively. Assuming 
the amplitude at is independent with target state Xk, we can rewrite L 2 (x) and c(z) as 

L z( x ) = L z( x )sl{a\d) (22) 

c{z) = c{z)c T a {a) (23) 

where L z (x) is the measurement location likelihood function and c(z) is the probability density of the 
false alarm spatial distribution in the image plane. We assume that the targets are within the surveillance 
region of sensor, the probability of detection for a given threshold x is then only dependent on d 

Pd( x ) = Pd(A) (24) 

Substituting Equations (5), (8) and (22-24) into Equation (20) we have the pseudo-likelihood 
function of AI-PHD as 

hM ) = l- P Ud) + P l{d)^ P uTnU [ t^ n U 1 ^ 

isz t A-V-c(z) + p D (d)-p a (a\d)-D kl[k _ x [_L z \ 

Equations (18) and (25) compose the recursion of AI-PHD filter. The probability of detection 
p T D (d) and amplitude likelihood ratio p T a (a\d) are replaced by p T D and p T a {a) respectively for the 
unknown SNR case. 

We can simplify the computation of p r a (a) by noting the fact that p T a (a) is calculated combined 
with p T D (d) in Equation (25). From Equations (8) and (9) we have 

p r D {d)-pl{a\d)= u (26) 
c a (a) 

Hence, instead of computing p T a (a) by Equation (13) we can compute the expression ■ p T a (a) 
directly using the method introduced in Section 2.2, i.e., 

• d 2 p x (a\v) 1 



c K a ) <?-c T a (a)-(d 2 -d 1 ) 



a ) _ I , 



Ol d 2 -O d y 
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(27) 
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where ®(x) = -j==j x exp -— dt is the standard normal distribution function which can be computed 



easily. Consequently, our approach incorporates the amplitude information into PHD filter with only a 
minor additional computational load. 

We show the consistency of our AI-PHD filter with the generic PHD filter. If the SNR of target is 
set as d = 0, from Equations (7-9) we have p T a {a\d) = \. This is the condition under which our AI-PHD 
filter degenerates to the generic PHD filter. 

3.3. Gaussian Mixture Implementation 

An analytic solution to the PHD filter can be found under linear assumptions on the system and 
observation equations with Gaussian process and observation noises as described in Equations (16) 
and (17) [20]. In this case, both the prediction and update equations of PHD are represented by a 
mixture of Gaussians where the means and covariances are updated with Kalman filter and the weights 
of the Gaussian components are found using the PHD filter equations. We use Gaussian Mixture 
implementation of our filter for its simplicity in calculation and convenience of target state extraction 
comparing to the Sequential Monte Carlo (SMC) method [21]. 

We assume that the survival probability is state independent, i.e., p sk {x') = p Sk and the detection 
probabilities are p T D (d) and p T D for the known SNR and the unknown SNR case respectively. The 
intensity of the target birth RFS are Gaussian mixture of the form 

y t (x) = £^N (*;««, P«) (28) 

i=l 

where J rk ,m^ k ,mf k ,P^ k ,i = \,---J^ k , are given model parameters that determine the shape of the 
birth intensity. 

We assume a uniform location distribution of clutter in the measurement space, so that the clutter 
location likelihood is not dependent on the state or the measurement. Hence the clutter location 
distribution is constant over the measurement space and equals to the reciprocal of the area of the 
image plane of optical sensor, i.e., c{z) = l/V. Next we give the prediction and update equations of 
Gaussian mixture implementation of our AI-PHD filter. 

Prediction: The posterior intensity at time k - 1 is a Gaussian mixture of the form 

A-i| i -iW = S4 ! 2iN (*;*««,/>«) (29) 

where J k -\ is the number of Gaussian terms with the weights w [ k \, means mf_ x and covariances P k %. 
Then the prediction intensity is still a Gaussian mixture as 

Jt -' 

*Vi(*) = fit(*) + />s,*5M-i N ( x ' m sM-v p s%k-i) (30) 

where the birth intensity y k (x) is given by Equation (28) and the means mf Mk _ x and covariances P^}^ 
are computed with the Kalman filter prediction. 
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Update: We rewrite the predicted intensity D k ^ k _ x {x) as a Gaussian mixture of the form 

Jk\k- \ 

i=i 



(31) 



Substituting the Equations (31), (20-25) into Equation (19), we obtain the intensity of our AI-PHD 
filter updated by measurements set Z k as the Gaussian mixture form 

D k]k (X) = [\-pl (d)] D k]k _, (x) + X X 4° (z)N (x; «$ (z), ) (32) 

where the updated means m^ k (z) = m ( k f k (z) and covariances P$ are calculated with the Kalman filter 
update. The updated weights o) k \z) in Equation (32) are computed as 

pUd)-pl{a\d)-a^_^ (z;z[ ! ,L„^;l-i) 



(33) 



where = .H^wi^ is the predicted measurement and =H k Pj^ c _ i Hl +R k is the innovation 
covariance. 

In the PHD update equation, Equation (32), and the weight update Equation (33), we replace the 
probability of detection p T D (d) and term p T D (d)- p T a {a\d) by p T D in Equation (14) and p T D -p T a {a) in 
Equation (27) when the SNR is unknown. 

4. Simulation 

In this Section, by setting up multi-target tracking simulation in the image plane of optical sensor, 
we examine the performance and computational complexity of our method for known and unknown 
SNR cases and benchmark them with generic PHD filter with different combinations of probability of 
false alarm and SNR value. 

4.1. Simulation Scene and OSPA Metric 



Consider a scenario with an unknown and time varying number of targets in clutter in the image 
region [-300, 300] x [2,000, 2,600] (pixel). Up to Nu = 6 targets are generated in this region with the 
random birth and dieing time instants. Figure 1 shows the true trajectories of each target. All targets in 
each simulation had the same mean SNR (this is not necessary by the algorithm but simplifies the 
presentation of results). 

Each target has survival probability ps,k = 0.99 and follows the linear Gaussian dynamics in 
Equation (16) with: 



k-\ 



I 2 AI 2 



*-i *i 

4 2 2 2 



I 2 A z / 2 



(34) 



where A = 1 s is the sampling period, and c v = 0.5(pixel/s ) is the standard deviation of process noise. 
The location measurement follows the observation model in Equation with H k =[I 2 0 2 ], R k =o- 2 £ I 2 , 
where c £ = \{pixel) is the standard deviation of measurement noise. 
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Figure 1. Target trajectories in the pixel plane with start/stop position as O/A. 
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The intensities of the birth RFS are Gaussian mixtures of the form: 



n (x) = f j o.m ( X \mf,p®) 

1=1 



(35) 



where mf 



is chosen around the mean initial states of i-th target, P^p = [50,50,15,15], i = l, - "6. To 
mitigate the exponential growth of mixture components, at each time step the number of Gaussian 
components is capped to a maximum of J max =100 components, whilst pruning is performed with a 
weight threshold of T= 10~ 5 , and merging is performed with a threshold of U= 2. 

We adopt the optimal subpattern assignment (OSPA) metric [22] for the purpose of multi-target 
performance evaluation since it jointly captures errors in the target state and target number estimates. 
Given two arbitrary finite sets X = {x 1 ,---x m } and Y = [y u - • • y n } , the OSPA is computed as follows: 

0 ,m = n = 0 



df{Xj) = 



V 



min V d 



(c \x i ,y m Y+c p {n-m) 



,m> n 



Up 



,m<n 



(36) 



where U n is the set of permutations on {\,-"n}, cf c \x,y) = min(c, d(x,y)), p is the order that penalizes 
error of individual element estimates, c is the cut-off parameter that penalizes error of cardinality 
estimate. We chose p = 2 and c = 30(pixel) in our simulation. Note that the chosen value of cut-off 
parameter c is significantly larger than the typical measurement noise but significantly smaller than the 
maximal distance between targets, thus maintaining a balance between the cardinality and localization 
components of the OSPA error [22]. 

4.2. Numerical Results 



4.2.1. Filtering Results for Multi-Target Tracking 



The effectiveness of our AI-PHD filter for multi-target tracking in image plane of optical sensor is 
verified through simulation. We assume a moderately cluttered scenario that the probabilities of false 
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alarm, p T F4 = 1 x 10~ 4 , which means the clutter density, X = 1 x 10~ 4 pixef 2 . The SNRs of all targets are 
set as J = 6 and the probability of detection p T D (d) ~ 0.99 (see Table 1). For the unknown SNR case, the 
SNR region is set as [2,10] and the probability of detection is replaced by p T D which can be computed by 
Equation (12). Other parameters for the filter are given as in Section 4.1. The true trajectories and filter 
estimates are shown in x and y coordinates of image plane versus time for AI-PHD filter with known 
and unknown target SNRs in Figures 2 and 3 respectively (denoted as easel and case2 accordingly). 

Figure 2. Filter estimates with known SNR. 
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Figure 3. Filter estimates with unknown SNR. 





300 n 


CD 
X 


200- 


te (pi 


100 


cc 
c 


o- 


coord i 


-100 
-200 


X 


-300^ 



■ True tracks AI-PHD filter estimates-case2 



Clutter 




Kbx|#fe 



*> — . 

_ 

1 



10 20 30 40 50 60 70 80 90 

time (s) 

(a) 



100 



Sensors 2012, 12 



2930 



Figure 3. Cont. 




From the estimates of AI-PHD filter shown in Figures 2 and 3, we see both for the known and 
unknown SNR case, the filters can eliminate dense clutter in the scenario with time -varying number of 
targets. All targets are detected immediately after birth and tracked accurately, which highlights the 
track initialization, maintenance and termination capabilities of our algorithm. The results also 
demonstrate superior performance of the algorithm in targets number and target states estimation. We 
evaluate this performance by Monte Carlo simulation next. 

4.2.2. Monte Carlo Results and Analysis 

Average OSPA and computation time cost per frame are used to evaluate the performance and 
computational complexity of our AI-PHD filter. To prove the improvements of the proposed method, 
we have benchmarked the results against a generic PHD filter which does not use the amplitude 
information. To make the assessment as fair as possible, the probability of detection p T D for filters 
without amplitude information was chosen to be the same as AI-PHD filter of the known SNR case, as 
in Table 1. 50 Monte Carlo runs were carried out for each combination of p T FA and d on computer 
(Intel quad core processor 2.66 GHz, 32-bit operating system, 4 Gbytes RAM) using Matlab (R2009b). 
Average OSPA for AI-PHD filter of known SNR and unknown SNR case and generic PHD filter are 
given in Table 2 where the results are divided by 7' accordingly. 



Table 2. Average OSPA (pixel) for different algorithms. 
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0.97/2.69/5.27 
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4 3.7190 


17.82/16.13/19.19 


7.37/6.44/11.01 
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1.11/3.41/7.82 
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16.41/13.53/20.21 
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2.79/4.14/14.89 


1.38/4.18/15.15 


0.95/3.99/14.43 


1 


x 10" 


3 3.0902 


14.82/11.83/21.42 


5.92/5.85/17.81 


2.26/5.54/17.37 


1.11/4.38/17.09 


1.00/4.03/16.94 



From Table 2 we see that for all combinations of p T FA and d given(corresponding to different 
probabilities of detection), our AI-PHD filter both with known and unknown SNR gives better 
performance than the generic PHD filter. This improvement is enhanced as p T FA or d increases. In the 
case ofp} A =1 x 10~ 3 and d = 8 where the method using the amplitude information works best, our 
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AI-PHD filter achieves 15.94 and 12.91 lower average OSPA (pixel) for known and unknown SNR, 
respectively. This improvement in performance is mainly due to two reasons: firstly, as d increases, the 
false alarm distribution will poorly represent the target counterpart and there is a big distinction in the 
target and false alarm distributions; secondly, as p T FA increases, having more measurements aids the 
method using the amplitude, since we discard less useful information. Table 2 also shows that the 
performance of generic PHD filter without amplitude decreases rapidly as p FA increases since there are 
more measurements from false alarms which by no means could be identified from those from targets. 
In contrast, we see no deterioration in the performance of AI-PHD filter in this case. For the known 
SNR case especially, the performance increases consistently as p T PA increases, which means our 
method works even better in a scenario with dense clutters. The comparison of computational 
complexity between AI-PHD filter and generic PHD filter without amplitude information is shown in 
average computation time per frame versus target SNR for different p T FA in Figure 4. Since we can 
achieve similar complexity for unknown SNR case with that of known SNR by computing 
Equation (27) with some fast algorithms, only the result for the known SNR case is given. 

Figure 4. Average computation time per frame for different algorithms with different 
SNRs d and probabilities of false alarm p FA . 
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We see that for different given p T FA and d, the results from two filters are close with the maximum 
difference being no more than 1 s. Figure 4(a) shows that in the scenario with low clutter density, 
AI-PHD filter has only a minor increase in average computation time over the generic PHD 
counterpart. Furthermore, AI-PHD filter performs an even low value in high clutter density scenario 
which is shown in Figure 4(b). In the case of p} A = 1 x 10~ 3 , d = 8 where this reduction is most 
obvious, the average computation time of AI-PHD filter is reduced by 53.7% over the generic PHD 
counterpart, which means the AI-PHD filter has even lower computational complexity than the PHD 
filter without amplitude information in scenarios with dense clutter and high SNR. The primary reason 
for this trend is that in these scenarios, the computation time cost is mainly decided by the multi-target 
state extraction step given the same number of targets and measurements. Incorporated with amplitude 
information, the update for the AI-PHD filter (see Equation (33)) gives heavier weights to the 
Gaussian items updated by the measurements from targets, thus updating the PHD with comparatively 
higher intensity near the real target positions and at the same time, suppressing the intensity of PHD 
near clutter positions (see Figure 5). Therefore, the updated Gaussian items can be prune and merged 
quickly and accurately. 

Figure 5. Intensity functions for generic PHD filter (left) and AI-PHD filter (right). 
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5. Conclusions 



In this paper, we have proposed a FISST-based method using signal amplitude information in a 
PHD filter for multi-target tracking applications in the image plane of optical sensors. We extend the 
measurement model to include the signal amplitude of observations and then incorporate this 
information into a PHD recursion in the form of an amplitude likelihood ratio. Based on the 
assumption that the amplitudes of the measurements from true targets are stronger than those from 
clutter, we show our method can significantly improve the performance over the one without 
amplitude information. Furthermore, simulation results also demonstrate that our method has much 
lower computational complexity in the scenario with high SNR and dense clutter, which makes sense 
for its practical implementation. Future work will involve improving the tracking performance for the 
targets with much lower SNR. 
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